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Abstract. There has been a growing evidence for the existence of magnetic fields in the 
extra-galactic regions, while the attempt to associate their origin with the inflationary epoch 
alone has been found extremely challenging. We therefore take into account the consistent 
post-inflationary evolution of the magnetic fields that are originated from vacuum fluctua¬ 
tions during inflation. In the model of our interest, the electromagnetic (EM) field is coupled 
to a pseudo-scalar inflaton cj) through the characteristic term 4>FF, breaking the conformal 
invariance. This interaction dynamically breaks the parity and enables a continuous produc¬ 
tion of only one of the polarization states of the EM field through tachyonic instability. The 
produced magnetic fields are thus helical. We find that the dominant contribution to the 
observed magnetic fields in this model comes from the modes that leave the horizon near the 
end of inflation, further enhanced by the tachyonic instability right after the end of inflation. 
The EM field is subsequently amplified by parametric resonance during the period of inflaton 
oscillation. Once the thermal plasma is formed (reheating), the produced helical magnetic 
fields undergo a turbulent process called inverse cascade, which shifts their peak correlation 
scales from smaller to larger scales. We consistently take all these effects into account within 
the regime where the perturbation of cf) is negligible and obtain Refr ~ 10~^®G, indicating 
the necessity of additional mechanisms to accommodate the observations. 
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1 Introduction 


Magnetic fields are observed on many different scales in our universe. In particular, it has 
been known that galaxies and galaxy clusters have their own magnetic fields with the typical 
strength, 10“®G [1-5]. Although their existence is observationally confirmed, the origin of 
these magnetic helds is still one of the open questions in cosmology. Recently, the gamma- 
ray observations from blazars suggest the existence of cosmological magnetic fields even in 
voids [6-13]. They set the lower bound on the effective strength of the void magnetic fields, 
given as [10, 14] 


Reff > 




( 1 . 1 ) 


where B is the present field amplitude with the correlation length A.^ The recent simulation 
indicates that if magnetic fields whose comoving strength is larger than 10“^®G exist prior 
to galaxy formation, they can be amplified to 10“®G as the disk of galaxy develops [15]. 
Therefore the reported magnetic fields in voids can seed those in galaxies and galaxy clusters. 
Hence “primordial magnetogenesis” can be a key ingredient of cosmology to explain the 
existence of magnetic fields in both galactic and extragalactic scales. 

Magnetogenesis during inflation has attracted attention, and many models have been 
proposed so far [16-30] (see also ref. [31]). In such scenarios, large-scale magnetic fields 
are expected to be generated during inflation via its causal production of fluctuations. Re¬ 
cently, however, the obstacles of inflationary magnetogenesis have been recognized, and it 
is understood to be quite difficult to generate magnetic fields only during inflation with 
the field strength satisfying the observational lower bound given by eq. (1.1) [32-37]. Since 
magnetic fields are substantially diluted after inflation due to the cosmic expansion, their 

^If one assumes that the suppression of the cascade emission is caused by the time delay of the cascade 
photons, one obtains a more conservative bound as BeS > [10]. 
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energy density needs to be extremely large during inflation, altering the inflation dynamics 
and/or inducing an additional curvature perturbation which exceeds the observed value in 
the cosmic microwave background (CMB) radiation [38]. Thus inflationary magnetogenesis 
faces three obstacles: (i) substantial dilution of the produced magnetic fields after infla¬ 
tion, (ii) too large electromagnetic energy spoiling inflation itself, and (iii) induced curvature 
perturbations being inconsistent with the CMB observation. 

In this paper, to overcome these three obstacles, we consider the axial coupling model in 
which a (pseudo-)scalar field is coupled to the electromagnetism through a term, 
where is the field strength tensor of the 17(1) gauge field, and is its dual. This 
coupling naturally arises under symmetries as the leading interaction of a pseudo-scalar and 
excites either of the left- or right-handed polarization modes of the gauge field. Helical 
magnetic fields, which are characterized by an imbalance between their two polarization 
states, have an interesting feature that part of its energy on small scales is transferred to 
that on large scales by non-linear effects of turbulent plasma with a high electric conductivity. 
This phenomenon is called the inverse easeade proeess. In practice, this process amplifies 
the comoving amplitudes of helical magnetic fields on large scales, and, as a result, their 
characteristic scale can grow. Therefore it can alleviate the substantial dilution of such 
magnetic fields on large scales in the cosmological context [39, 40]. Recently, the possible 
signatures of cosmological helical magnetic fields are reported in ref. [41, 42]. 

In the previous works on the axial coupling models, only the slow-roll regime of 4> 
was considered [18, 24, 26, 27]. To evaluate the generated magnetic fields, they applied an 
approximated analytic solution of the gauge field which is valid only during the slow roll 
regime without its backreaction onto the inflationary dynamics. However, we find that the 
generation of the electromagnetic fields is much more efficient at the end of the slow-roll 
regime and during the subsequent oscillation phase of (/>. Since the slow-roll approximation is 
no longer valid in these regimes, we numerically solve the equations of motion of these fields, 
consistently taking into account the backreaction from the gauge field to the dynamics of cf) 
and the background expansion. As a result, much stronger magnetic fields can be obtained 
than the ones in the preceding studies with the same setup [18, 24]. 

In our scenario, the electromagnetic fields are mainly generated right after the end 
of inflation, and hence we do not need to worry about the electromagnetic energy density 
spoiling inflation. Furthermore, since the peak scales of the produced magnetic fields are far 
smaller than the CMB scales, they are free from the CMB constraints. Therefore our scenario 
overcomes the three obstacles of inflationary magnetogenesis, and relatively strong magnetic 
fields can be generated. Unfortunately, however, in the case where (p is the inflaton in large- 
held models, the strength of the resultant magnetic helds does not reach the observational 
bound, eq. (1.1), by a few orders of magnitude. 

Two additional issues concern us in the scenario. One is the effect of electrically charged 
particles. We investigate the generation of magnetic helds by solving the coupled system of 
the background inhaton and the gauge held. However, if a sufficient amount of charged 
particles are present at the time of the magnetic held production, they may change the 
dynamics of the gauge held in a non-trivial way [43]. To avoid such cumbersome complexity, 
we require the inhaton decay rate is suppressed such that the charged particles are negligible 
during the production period. Nonetheless, since the reheating is delayed, the dilution of 
the produced magnetic held before the onset of the inverse cascade process is inevitably 
signihcant. The other issue is the non-linear effect of the axial coupling on the dynamics of 
the inhaton and of the gauge held. With this non-linear interaction taken into account, the 
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gauge field produced by the background inflaton inverse-decays to the inflaton’s fluctuations, 
which then in turn produces the gauge field. For the large value of the axial coupling constant, 
this successive non-linear process would change the dynamics of the inflaton and gauge field as 
shown in ref. [44]. Since fully handling the non-linear dynamics requires lattice simulations, 
we only consider the coupling constant small enough to ignore the effect of the inflaton 
perturbation on the gauge field dynamics. 

The rest of this paper is organized as follows. In sec. 2, we introduce our setup and 
review the approximated solution for the gauge field in the slow-roll regime. In sec. 3, we 
show the result of our numerical calculation and describe the dynamics of the model. In 
sec. 4, we take into account the inverse cascade process as the conservation of the helicity 
of the produced magnetic fields. Following the evolution after reheating, we exhibit the 
computation of the present values of the field strength and the correlation length for a given 
coupling constant in the model. In sec. 5, we discuss the effects of charged particles on 
reheating in various cases. In sec. 6, the amount of the inflaton perturbation generated by 
the inverse-decay of the gauge quanta is evaluated, and the strongest reliable value of the 
effective magnetic field strength in our frame work is presented. Finally we conclude in sec. 7. 

2 Model 

In this paper we consider the axial coupling model with the following Lagrangian: 

£ = - ^^<i>F,,F^\ ( 2 . 1 ) 

where (j) is the inflaton ,— dyA^ is the field strength tensor of the 17(1) gauge 
field A^, and F^^^ = Fp^/ {2y/—g) is its dual. Inspired by the axion case in which (j) is 
a pseudo-Nambu-Goldstone boson coupled to the axial vector current of charged fermions 
and the (pFF coupling arises from the chiral anomaly, we introduce the coupling constant 
as, a/f where a is a dimensionless constant and / corresponds to the axion decay constant. 
However, since this coupling is realized simply by symmetry consideration, we concentrate 
on phenomenology of the above Lagrangian without assuming any particular particle physics 
model. For simplicity, the potential is assumed to be quadratic form,^ 

V{(j)) = m = 1.9 X lO^^GeV. (2.2) 

In the flat FRW universe with the Einstein gravity, one can derive the equations of motion 
for the background inflaton (/>o(t) and the mode function of the gauge field A±{k,t) as [46] 

Mt) + ^HMt) + mVo(i) = j{E-B), (2.3) 

a k ■ 

T 7-4 

fa 

^ The recent Planck result [45] disfavors the quadratic potential of inflaton. In this work, however, we 
use it as an example to depict the consistent generation and evolution of the magnetic fields from the ifiFF 
term. Moreover, as we will show later, the dominant contribution to the magnetic fields in this model comes 
from the modes that cross the horizon near the end of inflation, and their amplification due to parametric 
resonance occurs during the early stage of inflaton oscillation. During this period, the quadratic potential is 
expected to be a good approximation of less simple forms, justifying its use. 


A±{k, t) = 0, 


(2.4) 


^2 

A±{k, t) F[A±{k, t) ^ 
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with 


(2.5) 


, „ 1 /" k d r, , ,9 , , ,91 

{E-B) = -^ j j^-a-[\A+\ - 1 ^- 1 ] 

where the overdot denotes the cosmic time derivative, and a is the scale factor. In this paper, 
we take the Coulomb gauge, = diAi = 0. We decompose and quantize the gauge field as 


Ai{t,x) = ^ 

A=± 


Sk 

(27r)3 


gifc-a=g(A)(^) 


a^^^Ax{k,t) + Al{k,t) 


( 2 . 6 ) 


where e\^\k) are the right/left-handed polarization vectors which satisfy €ijikje^^\k) = 
^ike[^\k), and are the creation/annihilation operators which satisfy the usual 

commutation relation, = {27r)^6{k + k')5^^. With the definition of the electric 

and magnetic fields in terms of the gauge field, Ei = —a~^dtAi and Bi = a~‘^eijkdjAk, one 
can derive eq. (2.5). 

Eq. (2.4) implies that either of two polarization modes A± has the tachyonic instability 
for k/a < a\4>o\/f, depending on the sign of (pQ. Thus the gauge fields are produced at 
the expense of the time kinetic energy of <po. The term on the right hand side in eq. (2.3) 
represents the backreaction from the gauge field to the inflaton and it slows down the motion 
of 4>q. It should be noted that the perturbation of the inflaton, 5(p{t,x), is ignored in the 
above set of equations. We discuss the validity of this assumption in sec. 6. 

Before exploring our scenario with the full numerical calculation, it is instructive to 
review the gauge field production in the slow-roll regime of p. For this purpose, it is useful 
to introduce a new parameter ^ as 


_ apo 
^ VH' 


(2.7) 


where H is the Hubble parameter. With the conformal time r], eq. (2.4) is rewritten as 


92 + ± 2k- 

I 


= 0 - 


( 2 . 8 ) 


When is constant which is a good approximation in the slow-roll regime, the exact solution 
of this equation is available. If pQ < 0, then ^ < 0 and the (—) mode has the tachyonic 
instability. With the adiabatic initial condition, 

A±{\kr)\ > 1) = (2.9) 

y/2k 

in the sub-horizon limit, one finds the analytic solution and the asymptotic expression for 
the growing (—) mode are given by 


A-{k,r]) 


1 


Go (1^1, -krj) + iFo (|^1, -k-q) 




( 2 . 10 ) 


up to an irrelevant constant phase, where Fl{h,z) and Gl{h,z) are regular and irregular 
Coulomb wave functions, respectively. Therefore the growing mode is amplified by a factor 
2'k\P\ due to the tachyonic instability. Since the produced gauge field induces non- 
gaussian perturbations of the inflaton, 54>, the value of ^ at the horizon-crossing of the 
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Figure 1. The time evolution of oc ^ is shown. The horizontal axis denotes the e-folding 

number N, and inflation ends {ch = —H/H^ = 1) at = 0. The initial conditions are <f)Q = •\/282Mpi 
and 00 = — y^2/3Mpi without any gauge fields. The blue dashed line shows the case with a = 0, while 
the red line shows the case with aMp\/f = 8. In both cases, (t>o/H (or increases during inflation 
and reaches its maximum value after the end of inflation. Therefore the most efficient production of 
the gauge fields takes place after inflation. 

CMB mode (A:* = 0.002Mpc“^) is constrained as |^*1 < 2.37 [47]. At smaller scales, non¬ 
detection of primordial black holes (PBHs) can potentially push the constraint further down 
[48]. While the constraint from PBHs inherits uncertainty, the forthcoming generations of 
terrestrial gravitational-wave (GW) detectors can directly probe the GWs sourced by the 
produced gauge field. It has been demonstrated that, for a quadratic potential of inflaton 
0, the second-generation detectors (Advanced LIGO, VIRGO, and KAGRA) will be able to 
probe |.^*| > 2.2 (at 95% GL, based on flat priors on and the third-generation experiments 
can improve the sensitivity to > 1.9 [49, 50]. 

Since the gauge field production is characterized by |^|, it is expected that the generation 
of the gauge fields is the most efficient when |,^| reaches its maximum value. Indeed, ^ oc 
0 o/-H ~ —2Mpj/0o in the slow-roll regime, and its absolute value increases during inflation. 
In fig. 1, we plot the evolution of (po/H oc One can see that |^| increases during inflation 
and reaches its maximum value when 0o crosses its origin at V ~ 0.5 after inflation. It 
then starts to oscillate with a large amplitude in the oscillation phase. Thus, in order to 
evaluate the produced magnetic field, one has to investigate the dynamics after inflation, the 
period during which the analytic solution is no longer valid. It seems that this critical point 
is overlooked in the previous works [18, 24, 26, 27]. Note that the amplitude of ^ slightly 
decreases and the oscillation phase changes for V > 1.5 in the case with aMpi/f = 8, because 
of the backreaction from the gauge field (compare the two lines in fig. 1). 

3 Numerical Result 

By a numerical calculation, we simultaneously solve eqs. (2.3), (2.4) and the Friedmann 
equation. 




+ Pern (a) 


(3.1) 




















where M-p\ is the reduced Planck mass, and the energy density of electromagnetic helds pem 
is given in terms of the mode function of the gauge field A±{k,t) as 


Pern — P 


(-) 


+ P 


(+) 
em ! 



k^\A±\^ 

27r2a2 ’ 


Pern ~ 2 k ^ ^ ^ ' B 

,(±) ^ k^\A±\'^ 

B 27r2a4 • 


(3.2) 

(3.3) 


However, since the gauge field is solved in Fourier space, and {E ■ B) and pem contain 
the momentum integrations, we have to solve eq. (2.4) for many modes of A±{kA) with 
different wave numbers simultaneously. In our numerical calculation, we start to solve for 
a mode function A±{k,t) with the initial condition eq. (2.9), when its wave number equals 
to a rapidly growing function, lOOfeins(t), where fcins = aacj)/f = 2^aH is the wave number 
for which the (—) mode becomes unstable, and continue to solve the mode until the end of 
calculation. In other words, we solve for all modes which satisfy 100A:ins(ii) < k < 100k^^^{t), 
where ti is the onset of the calculation, N{ti) ~ —70, and is the maximum value of 

kins on the interval [ti, t]. The grid size of the wave number is set as din /c = 0.01. To eliminate 
the contribution from the modes in the vacuum state, we do not integrate all the modes being 
solved in eqs. (2.5) and (3.2) but only the modes satisfying k < hnsit) at least once during 
the numerical evolution. For conhrmation, we check another criterion. We integrate only 
the modes with a sufficient amplitude, |\/^.A±| > 20, and compare the result with that 
of the former criterion, confirming no recognizable change. The time step of the numerical 
calculation is set as 6t = 10“^a(f)/fcins(t), which is always much shorter than the inflaton 
oscillation time scale, m~^. We check the total energy conservation, p = —3H{p + p), is 
satisfied in sub-percent accuracy. 

In this section, we discuss only the case of aMpi/f = 8, because it is close to the 
maximum value of the coupling constant aMpi/f for which our treatment is justihed. For 
larger aMpi/f, we need to take into account the inflaton perturbation induced by the gauge 
held. However, it requires lattice simulations as the re-scattering process between the gauge 
held A± and 64> is fully non-linear. Thus in this paper, we only consider the case where 
the coupling is not too strong and the inhaton perturbation is negligible, 54> ^ (po- We will 
discuss this point with more detail in sec 6. 

Let us now show the result of our numerical calculation. In hg. 2, we plot the energy 
density of the inhaton and the two polarization mode of the gauge held in the left panel. 
The horizontal axis N denotes the e-fold number, and inhation ends at = 0. One can 
clearly see that the electromagnetic helds are mainly generated after inhation and the energy 
density pem = Pem^ + Pein^ reaches its maximum value at 2 e-folds after the end of inhation. 
Note that pem is always dominated by the contribution from the (—) mode pim . Until 
N ~ 0.75, cp is negative (see hg. 1), and only the (—) mode of the gauge held continues to 
grow due to the tachyonic instability. However, after the sign of (p changes, the (—) mode 
transiently decreases due to acquiring a positive mass, while the {+) mode begins to increase 
by the tachyonic instability in turn. Subsequently, the inhaton oscillation develops, and the 
parametric resonance becomes efficient. Although not only A- but also A+ is amplihed 
in this regime, A- is far larger than A+ due to the initial hierarchy implemented by the 
tachyonic instability, and thus a helical magnetic held is generated. As the amplitude of the 
inhaton oscillation decreases, the efficient resonance eventually terminates. After that, since 
the inhaton and the gauge held are now decoupled, they behave as the usual matter and 
radiation components, respectively. 
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Figure 2. The energy density of the inflaton (blue line), the (—) mode of the gauge Held p^~'^ (red 
line) and the (+) mode of the gauge field p^~^'> (yellow line) are shown in the left panel. The potential 
force V' = dVjdcj) = (blue line) and the backreaction term a{E ■ B)/f (red line) in eq. (2.3) 
are compared in the right panel. We take aMpi/f = 8. The tachyonic instability and the subsequent 
parametric resonance amplify the (—) mode while the (+) mode is not significantly produced. The 
energy densities eventually decay as pem cx; a~^ and oc a~^ after the coupling between (j) and A± 
becomes ineffective at ~ 2.3. 

yjAA. yJYt A, 




Figure 3. The time evolution of the spectra of the mode function A±{k,t)- The horizontal axis 
denotes the physical wave number at the end of the calculation, Nf Ri 2.9. Colored lines represent 
the spectra at different times with time steps of 1 e-fold from bottom to top, and the top red lines 
denote the spectra at Nf. A- (left panel) acquires the double amplifications both from the tachyonic 
instability and the parametric resonance, and the peak mode reaches amplitudes of 0(10®). A+ (right 
panel) is slightly damped during inflation due to the positive effective mass, while the parametric 
resonance amplifies the peak mode by O(IO^). 


In the right panel of fig. 2, the potential force of the inflaton V' = dVidcj) = m?cj) 
and the backreaction from the gauge field to the inflaton a{E • B)/f are shown. Although 
the backreaction term is always subdominant in this case of aMp\/f = 8, it becomes non- 
negligible and the behavior of (j) is significantly affected for a larger coupling constant. 

In fig. 3, we show the spectra of the mode function A±{k,t)- Although A+ does not 
have the tachyonic instability during inflation, the oscillating inflaton amplifies A+ to at 
most O(IO^) through the parametric resonance as the oscillating behavior of the spectrum 
indicates. The mode A- acquires the double amplification both from the tachyonic instability 
and the parametric resonance, and the peak mode at k/af ~ 0.1m undergoes the 0(10®) 
amplification. It corresponds to |i^| ~ 4.5 in the constant ^ case, while ^ at the horizon 
crossing of the CMB scale is |^| ~ aMpl/2^/ANcMB'^^f ~ 0.3 in this setup. It should be 


- 7 - 














































































^phy ^phy ] 




Figure 4. The time evolution of the physical intensity (left panel) and the physical correlation 
length (right panel) of the magnetic field are shown. The e-fold number N is normalized such that 
inflation ends at = 0. After the decoupling (iV > 2), they behaves as oc a“^, and Aphy oc a 
(red dashed lines). However, this behavior will be changed when charged particles appear and begin 
to interact with the electromagnetic helds. 


noted that since A- is much larger than A+ at the end, an almost maximal helical magnetic 
field is eventually generated (see eq. (4.1)). 

Following refs. [27, 40], we define the physical intensity and the physical correlation 
length of the magnetic field as 



AphyW_a + 


(3.4) 


In fig. 4, we plot them as a function of the e-fold number N for aMp\/f = 8. After the 
coupling between the inflaton and the gauge field becomes ineffective {N > 2), the electro¬ 
magnetic fields effectively satisfy the free wave equation, [5^ -|- k‘^]A±{k,r]) = 0, and thus 
|,A±| become constant. Then we can fit the numerical result as (see the red dashed lines in 
fig. 4.) 

~ 1.3 X 10^ {N>2). (3.5) 

By using the fitted behavior of the Hubble parameter, H ~ 0.36me“^^/^, and the unit 
conversions IG = 6.8 x lO^^'^GeV^ and lGeV“^ = 6.4 x 10“^®Mpc, one can rewrite eq. (3.5) 
as 

8/3 . s -2/3 

G^, Aphy ^ 4.5 X lO-'^^ _ Mpc, [H < IQ-^m). (3.6) 


B: 


phy 


5.4 X 10^°^ 


- 

m J 


This is the physical magnetic field amplitude with the associated correlation length after all 
the production processes complete, and this behavior continues to hold until charged particles 
are generated by the inflaton decay. Thereafter the produced helical magnetic field undergoes 
the inverse cascade process. 


4 Inverse Cascade Process 

As the inflaton decays and particles in the standard model appear, the electric conductivity Cc 
increases. The conductivity gives a friction term in the equation of A± as [d^ -|- {H + ac)dt + 
k"^/a^]A± = 0. Since, after the reheating completes, the conductivity of the thermalized 
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plasma exceeds the Hubble parameter, ac^ H, one can find that on large scales k <C ay/acH 
the gauge mode functions freeze, A± ~ const. [25]. Thus the corresponding electric fields 
vanish, and the static magnetic fields which merely dilute due to the cosmic expansion remain 
on such large scales. On small scales, the evolution of magnetic fields is different due to the 
nonlinearity of the fluid dynamics [39, 40]. Turbulence in a plasma induced by the Lorentz 
force highly develops on small scales. The turbulence transfers the energy of the magnetic 
fields from large scales to small scales due to a non-linear interaction between the plasma fluid 
and the magnetic helds. This is called the direct cascade process. The resultant spectrum 
of the magnetic field is red-tilted on small scales where the time scale of the turbulence is 
smaller than the cosmological time scale. 

However, when a magnetic field is helical, not all of its energy is transferred to small 
scales. According to Maxwell equations, the magnetic helicity should be conserved in the 
magnetized plasma with a high conductivity. Therefore, in order to maintain the helicity, 
a part of the energy is transferred to a slightly larger scales where the turbulence does not 
develop well. This is called the inverse cascade process. 

The magnetic helicity % is defined as "H = jy d^xCijkAidjAk where the spacial integral 
is taken over a finite comoving volume V on which magnetic fields have no component normal 
to the boundary, i.e. B ■ n = 0. Then T-L is invariant under the gauge transformation (see 
e.g. [51]). With eq. (2.6), the expectation value of the helicity density f) which is defined as 
f) = 7i/V is given by 

{i]) = ^ J^d^xcijkiAidjAk) = J (4.1) 

Thus the helicity density represents the difference between the left-handed and the right- 
handed polarization modes. When a magnetic field is perfectly helical (i.e. one component 
of polarization modes is non-zero while the other components are zero), the helicity conser¬ 
vation, Ti = const, or equivalently 1) = const., provides the constraint between the strength 
and the correlation length of the helical magnetic field as 

(f)) oc a^AphyHphy = const. (4.2) 

It is clear that this equation holds for freezed magnetic fields because their physical scale and 
strength simply evolved as Aphy(t) oc a and Bp^y{t) oc a“^. However this equation is valid 
even in the turbulent regime where the magnetic field evolution is affected by the turbulence 
in a plasma and Aphy grows faster than a. The scale and strength of magnetic fields in the 
turbulent regime evolve with satisfying this equation. 

The study about the evolution of magnetic fields with non-zero helicity after inflation 
requires a MHD numerical simulation. However the helicity conservation enable us to directly 
compare helical magnetic fields produced in early universe with the observation bound by 
using eq. (4.2), because the constraint from blazar observations, eq. (1.1), is proportional to 
AphyHpj^y oc (1)) for magnetic fields with A < 1 Mpc. The scale factor at the completion of 
reheating is given by 

a,eh~ 5.5 X 10-32 = (4.3) 

where T^ is the decay rate of the inflaton, and we set the scale factor at present as unity. 
Therefore, when the helicity conservation is valid, the effective strength of the helical magnetic 
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(A < IMpc), (4.4) 


field at present can be written as 



where treh denotes the time at the reheating completion. Recall Reff is defined in eq. (1.1). 
Snbstituting eq. (3.6) with H = we finally obtain the produced Reff at the present epoch 
in the case of aMpi/f = 8 as 

B., = 6.6x10-»g(^)''‘, (^^=8). (4.5) 

In the next section, we discuss the allowed values of Furthermore, in sec. 6, we explore 
the larger values of the coupling constant a/ f. There exists the parameter space where the 
stronger Reff is obtained. 

Although, with the helicity conservation, we can estimate Rgg at present from the 
strength and scale of a helical magnetic field at the reheating epoch, Rphy(t) and Aphy(t) are 
degenerate in Reff. Before closing this section, we explore the way to resolve the degeneracy 
between Rphy(t) and Aphy(t). 

If the initial spectrum of the helical magnetic helds is blue-tilted, as in the present case, 
its characteristic length scale, at which most of the magnetic field energy resides, is shifted 
to the scale comparable to that of the turbulence, XT{t). This is due to the direct cascade 
process, in which the magnetic field energy on the scales smaller than Ar(t) is significantly 
transferred by the turbulence to farther smaller scales and dissipates away. 

The turbulence scale At is given by the plasma fluid velocity v and the cosmic time 
t as Xrit) ~ vt. When the turbulence fully develops, the equipartition state between the 
fluid and the magnetic field is achieved. In that case, v reaches the Alfven velocity, v^it) ~ 
Bphy/y/p, where p is the energy density of the plasma. Thus one finds Aphy ~ tRphy/y^. This 
equation holds until the recombination when most of the charged particles disappears and the 
turbulence terminates (see also refs. [39, 40] for a more detailed discussion). After that, Rphy 
and Aphy evolve as usual, only due to the cosmological expansion, Rphy oc a~^,Aphy oc o. ^ 
Thus, Rphy/Aphy at present is given by 


R, 


phy 


\phy 


(R 


■sj p(Rec) ^3 


P{trec) 

Mpi 


„3 -'CMB-^rec 

O, 

'^rec 


Mpi 


(4.6) 


where the subscript “rec” denotes the recombination and Tcmb ~ 3K. Therefore, we ob¬ 
tain [39, 40] 

Bpk,(t„ow) ~ 10-»G (hi^) . (4,7) 

Combining this equation and the helicity conservation, one can solve for each of Rphy and 
Aphy at present as 


Bpl,y(t..w) ~ 10-*A-G, Aphy(t.„„) ~ ATMpc, A- = V''” 

(4.8) 

^ Since residual charged particles remain after the recombination, the MHD interaction is effective between 
the magnetic fields and the plasma fluid. However, since the Alfven velocity is proportional to in the 

matter dominated plasma, tv a evolves as cx: a. Therefore, the MHD interaction in the matter dominated era 
does not induce the further evolution of Aphy/a. 
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Substituting eq. (3.6) with H = F^, we find 




3.5 X 10“® 


V 

lOSGeVy ’ 


/ aMpi 

v“r 



(4.9) 


The MHD numerical simulations in ref. [39, 53, 54] show that, due to the inverse cascade, 
the spectral index of the magnetic field energy spectrum, is ~ 5 on the scales larger than 
the peak scale. Therefore, the magnetic field strength at the scale L {L > Aphy) is roughly 
evaluated as Bl i?phy(Aphy/L)5/2. 


5 Reheating 

Provided that reheating proceeds instantaneously and charged particles appear at once, one 
can substitute eq. (3.6) with H = into eq. (4.4) independently on the value of T,^, (if 
A < IMpc and Treh ^ IO^^ttt,). However, charged particles can be produced before the 
completion of reheating, R = T^, since the inflaton decay is active even for R > T^. For 
example, by assuming the inflaton decay rate T^ is constant and the energy density of the 
decay products /Orad(i) satisfies, prad + dRprad = one finds prad = 2r0P(^/5R in the 

inflaton oscillating phase [52]. Thus the decay products occupy roughly V^/H of the total 
energy density. If the charged particles in prad affected the dynamics of A± before the 
decoupling between the inflaton and the gauge held, we would have to reconsider the result 
in sec. 3, which does not take into account the effect of the charged particles.'^ Therefore in 
this section, we discuss the maximum value of T^ for which our numerical result in sec. 3 is 
verihed. 


5.1 Case 1 : the inflaton decays into charged particles 


First, we consider the case where the inflaton mainly decays into charged particles with a 
constant decay rate F^. The produced charged particles alter the dynamics of the gauge 
held if the interaction between them is signihcant. Since the cross section of the interaction 
is roughly Uint — /]? where ol is the hne structure constant, and p denotes a typical 

momentum of the charged particle whose mass is negligible compared to the inhaton mass 
m, the interaction rate per a Hubble rate can be estimated as 

hint _ WcUintU ^ ''' f Y f ^ ^ ^ 

H ~ H ““ m3 ^VO-Oly VlO^GeVy V2 X 1033GeV 



where ric — Prad/w- is the number density of the charged particles whose momentum is p ~ m, 
and we have used prad — ~ Therefore, for F^ <C lO^GeV, the charged 

particles do not signihcantly interact with the gauge held within a Hubble time during the 
inflaton oscillating phase. Note that the electromagnetic flelds are generated for an interval 
of 2 e-folds after the end of inflation (see fig. 2). Thus F^ does not have to be suppressed by 
many orders of magnitude than lO^GeV. 

Eq. (5.1) is also used as a condition of thermalization, if a' is replaced by the hne 
structure constant of the coupling by which the charged particles reach thermal equilibrium. 
If the charged particles thermalize, the conductivity would be roughly given by [55] 

1/4 / N 1/4 


CTr. 


lO^T 


/ r \ ^ ^ / M \ ^ ^ 

lOVrid ~ 5 X lO^SQgy j ( iQiiGev J ’ (thermahzed). (5.2) 


'^In ref. [43], the authors study how the growing electric conductivity suppresses the generation of magnetic 
fields by a parametric resonance, while they assume different couplings, and RA^, from ours. 
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Hence cJc would be much higher than the Hubble at the decoupling between the inflaton and 
the gauge field, Hdec ~ lO^^GeV, and the generation of the helical magnetic fields would be 
highly suppressed. However, this estimation is invalid unless the interaction which brings 
the decay products into thermal equilibrium is effective. Tint > H. Thus, eq. (5.1) gives the 
relevant constraint on as 

r^«10'Gev(^) (5,3) 

Hence we substitute = lO^GeV in eq. (4.5). 


5.2 Case 2 : the inflaton decays into non-charged particles 

Second, we consider the case where the inflaton mainly decays into non-charged particles 
with a constant T^, while the non-charged particles can produce charged particles through 
an interaction with a cross section, (Tint — In this case, since the energy density of 

the charged particles is additionally suppressed by factor of fint/i? ^ m^, the 

interaction rate between the charged particles and the gauge held is given by 

Hint _ ncTintU ^ / a' ^ f ^ Y f ^ 

H ~ H “V"" m3 y VO-Oly VlO^GeVy V2 x 10i3GeV 

(5.4) 

The upper bound on Tj^^ is relaxed, since its dependence is squared. 

Provided that the inflaton (j) is an axion, namely as in the natural inflation case [56], it 
mainly decays into gauge fields and the decay rate is written as [47] 



4>^AA 


9 

am 


647r/^ 


m 






(5.5) 


where we consider the decay into the 17(1) gauge field for simplicity, while the axion may 
decay into other non-Abelian fields. Decays into fermions are further suppressed due to the 
helicity by a factor with a fermion mass < m [47]. Since the 17(1) gauge held 

does not have a self-coupling, eq. (5.4) is applicable to this case. Then, charged particles 
are unlikely to affect the dynamics of the 17(1) gauge held calculated in sec. 3. However, 
with such a small decay rate, the inflaton oscillation phase lasts for a long period and the 
produced magnetic fields substantially dilute as i?phy(t) oc a“^, until the reheating completes 
and the inverse cascade process becomes effective. 


5.3 Case 3 : the inflaton decays via Ynkawa conpling 

Finally let us discuss the case in which the inflaton mainly decays through a Yukawa coupling, 
where y is the coupling constant, and 'tp is a fermion. Although the perturbative 
decay rate of this coupling is Ty = y'^mjSTr, the actual decay process is non-trivial. First, 
the inflaton field value gives a large effective mass to the fermion, = yep, and the fermion 
is heavier than the inflaton until (p sufficiently damps. Second, during the cp oscillation, the 
fermion can be produced when the inflaton passes through its origin, (/> ~ 0, by the so-called 
fermionic preheating [57, 58]. Although the non-perturbative fermion production is studied 
in the literature, its backreaction to the inflaton is non-trivial in this process. Therefore it is 
beyond the scope of this paper to precisely calculate the maximum F ^ in the case of Yukawa 
coupling. 
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However, as a trial, let us estimate the decay rate by ignoring the fermion preheating. 
Provided the inflaton does not decay into the fermion nntil its oscillation amplitude (pose 
decreases to m/2y (namely = ml2), an optimized Ynkawa coupling y is obtained as 

where we have used H ~ mcposc /Note if y is larger than this valne, reheating 
is delayed because the fermion effective mass becomes larger and prevents the decay. The 
conpling constant in eq. (5.6) corresponds to the decay rate P^^ « lO^GeV. Therefore it would 
be interesting to further study the Yukawa coupling case, since it may realize a high reheating 
temperature avoiding charged particles from snppressing the magnetic held production. 


6 The inflaton pertnrbation 5(j) 


In sec. 3, we foens on the case of aMpi/f = 8. One expects if the coupling constant is 
larger, stronger magnetic helds are produced. However, at the same time, one should care 
about the consistency of the assumption, 6(j) (po, which we made in the eqnation of motion 
for the gauge held, eq. (2.4). If the coupling between the inhaton and the gauge held is 
too strong, the perturbation of the inhaton 6(p is signihcantly prodneed by the gange held 
and 6(p becomes non-negligible compared to (pQ. Then 5(p potentially alters the prodnetion 
of the gauge held, while it is difficult to take into account the non-linear coupling between 
6(p and A± without performing fnll lattice simnlations. To hnd the allowed maximum value 
of the coupling constant where our assumption is valid, we evaluate the variance of the 
inhaton perturbation {5(p‘^) by ignoring the backreaction from 6(p to the gauge held. That 
“maximum” coupling constant leads to the strongest magnetic held which is obtained under 
our assumption. 

The equation of motion for S(p{r], x) is given by 


d^-df- — + a^ 


2 , ^rn^(po4>o 

m + 2 - o - 


+ 


3 + 2 ^ ^ 


. a 


{ab(p) = a:^-[E-B-(E-B)], 


( 6 . 1 ) 

where the gravitational coupling to the backgronnd inhaton (pQ is also taken into account.^ 
The right hand side of eq. (6.1) represents the contribution of the gauge held through the 
coupling. The solution for the Fourier mode of 5<p{r], x) coming from the source can be found 
by using the Green function as 


a{rf)6(p{k,rf) =2 J dTlm[Q*{k,r])Q{k,T)] Jem{k,T) 


( 6 . 2 ) 


where Q{k,r]) is the homogeneous solution of eq. (6.1), and Jem is the Fourier-transformed 
source term, 

Jem{k,r]) = J ■^J^Ei{p,r])Bi{k - p,r]). (6.3) 

®Since these gravitational coupling terms cause the metric preheating of 5if> [59-61] and the oscillation 
amplitude of Scj> does not decrease while (po damps as after the decoupling. Thus Sep may become non- 

negligible eventually and contribute to magnetogenesis if the metric preheating lasts sufficiently long. We will 
come back to this issue in a future work. 
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Figure 5. The backreaction measure {5(f?)/p^j, at the time when Pem!P 4 > reaches its maxi¬ 
mum value (left panel), and the effective magnetic field strength B^s with = 10®GeV (right 
panel) are shown as functions of the coupling constant. Since p^/rn^ represents the cycle av¬ 
erage of (poit), p^ should be smaller than unity to ignore b(p in eq. (2.4). This con¬ 

sistency condition (aMpijf < 8.4) restricts the reliable prediction for BeS in our framework as 
BeS ^ 10“^®G (r^/lO^GeV)^^^. The fitted lines (red dotted) in the left and right panels are given in 
eqs. (6.7) and (6.8), respectively. 


After a few lines of algebra, we obtain the power spectrum of the 5(p induced by the gauge 
field as 


'Ps^jyik,?]) = 


a 




27:"^ pa?{ri) 


X 


E 

X,(t=zL 
2 


d^p 

(27r)^ 


1 — Xap ■ k — p 


P \^Xa{r,k]p,\k - p\)\ + Xap\k-p\Ixa{T,k-,p,\k - p\)I*^{T,k-,\k - p\,p) , (6.4) 


with 

Ixa{v,k-,p,q) = J -^lm[Q{k,p)Q*{k,T)]Ax{p,T)A'^{q,T). (6.5) 

The variance of the inflaton perturbation is given by {6(p‘^{r])) = f dkVs(j){k,p)/k. We numer¬ 
ically evaluate by using the solutions of cj)o{t) and A±{k,t) which are also numerically 

obtained in sec. 3. We compare {dcp'^) with which gives the cycle average of 0 q. 

In the left panel of fig. 5, we show rin?{5(lP‘) j at the time when the ratio between the 
energy density of the gauge field and that of the background inflaton, namely p^m/Pcft-, reaches 
its maximum value for various values of the coupling constant, aMpi/f = 7, 7.5, 8,8.3, 8.5,9 
and 10. For example, in the case of aMpi/f = 8, {6cj)‘^) plotted in fig. 5 is evaluated at 
~ 2. For aMpi/f > 8.4, the inflaton perturbation 6(j) becomes larger than the background 
value (pQ before the production of the magnetic fields terminates. Therefore our treatment is 
justified only for aMpi/f < 8.4. For aMp\/f = 8.3, we find m?‘{5(tP‘) jp^ = 0.2 when Pem/Pep 
reaches its maximum and then we obtain 

= (^ = 8 . 3 ). (6.6) 

Therefore this is the strongest magnetic field in our treatment. In the right panel of fig. 5, 
Bgf[ at present is shown as a function of the coupling constant, aMpi/f. One can clearly see 
that the resultant magnetic field becomes stronger as the coupling is larger. 
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We find that rin?{5(lP ‘)/and Sefr are well approximated by the following function (the 
red dotted lines in fig. 5): 

= 4 X 10-^» exp[7aMp,//], 

= 8 X 10"^'G exp [2aA/pi//] < gj (g gj 

For aMpi/f > 9, where our computation is invalid as 64> > 4’o, pem becomes much larger 
than p 0 in our numerical calculation, not only because the amplification from the tachyonic 
instability gets stronger but also because the backreaction from the gauge field suppresses the 
amplitude of the inflaton. Since pem P 4 >-, the universe is radiation dominated for a while 
after the decoupling, and the dilution of the produced magnetic fields is less signihcant. 
Accordingly, i?eff becomes stronger as the coupling is larger, although the almost all energy 
of the inflaton is transferred into the gauge held for aMpi/f > 9. However, such a large 
coupling constant invalidates the assumptions in our calculation, and the result is likely to 
be modihed, at least quantitatively, if 6(j) is properly taken into account. ® 

The result, eq. ( 6 . 6 ), for aMpi/f = 8.3 and F^ = lO^GeV corresponds to the magnetic 
helds with the strength Hphy(tnow) ~ 5 x and the correlation length Aphy(tnow) ~ 

0.05pc. Thus, in spite of the signihcant growth due to the inverse cascade, the correlation 
length of the magnetic held is still very small. It is because we have considered the large- 
held model of inhation in which the energy scale around the end of inhation is very high 
and the magnetic helds are generated on a very small scale. Note that the mass scale 
m = 1.9 X lO^^GeV corresponds to the length scale m~^ = 3 x 10“^^Mpc. 

In order to produce magnetic helds with a larger correlation length, one can consider 
a lower-energy inhation model or a spectator held with the axial coupling which begins to 
oscillate at a lower-energy scale after inhation. Although such a model may be less simple and 
contain more parameters than our setup in this paper, it is expected that one can calculate 
resultant magnetic helds in a similar manner. We intend to investigate this possibility in the 
future work. 

A spectator held with a generalized axial coupling in a low energy inhation has been 
considered in refs. [26, 27]. In these works, the authors assumed that ^ is constant until 
instantaneous reheating without solving the dynamics of the spectator held. In this paper, 
however, we have shown that ^ can vary in a non-trivial way and it leads to an efficient 
amplihcation of magnetic helds. We have also discussed the effects of charged particles and 
the perturbation of the axial coupled scalar held. It is also worth reconsidering these models 
taking them into account. 

® In the recent paper [44], the coupled dynamics of the inflaton and gauge field is studied by lattice 
simulations for larger coupling cases. The simulations with analytically approximated initial conditions of 
inflaton and gauge field at the end of inflation show that the helical asymmetry of the gauge field would 
be erased by the re-scattering process for a/f > 9M~^. However, the homogenous mode of the inflaton 
and the configuration of the gauge field deviates from the analytical form for the larger axial coupling as 
mentioned in [44], which has the possibility to change the eventual helical asymmetry. Thus, in order to 
evaluate the abundance of the magnetic field for larger coupling at least around q// ~ 9M~^, one needs to 
execute the non-linear simulations such as lattice simulations to take into account the re-scattering process 
with appropriate initial conditions. 
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7 Conclusion 


In this paper, we investigate the generation of magnetic fields in the axial coupling model 
where the inflaton and the gauge field are coupled through the term, C\at = 

Although only the slow-roll regime of the inflaton is considered in the previous works, we point 
out, for the first time in the context of magnetogenesis, that the most efficient production of 
the electromagnetic fields takes place after the end of inflation. This property of the model is 
quite interesting because it is known that magnetogenesis only during inflation suffers from 
several problems and thus the production of magnetic field with the sufficient strength to 
explain the blazar observations is difficult. 

We have numerically solved the coupled equations of the background inflaton and the 
gauge field by taking into account the backreaction from the produced gauge field to the 
inflaton and the background dynamics. Since the axial coupling spontaneously breaks the 
parity symmetry, one of the two polarizations of the gauge field is mainly produced. In 
our case where (p is initially negative, only the (—) mode, A-{k,t), acquires the double 
amplification from the tachyonic instability and the parametric resonance. For the opposite 
signature of p, only the (+) mode enhances. As a result, the almost completely helical 
magnetic field is produced. After turbulence of the magnetized plasma develops, the helical 
magnetic fields undergo the inverse cascade process, and their comoving amplitude grows. 
In virtue of the helicity conservation, we can easily calculate the effective magnetic strength 
Bgf[ at present from the physical intensity and correlation length of the magnetic fields at 
reheating. 

In this paper, we make two assumptions to simplify the calculation. We assume that 
electrically charged particles produced by the inflaton decay do not significantly interact with 
the gauge field until the production of the magnetic fields terminates, and that the coupling 
constant of the axial coupling is small enough to ignore the effect of the inflaton perturbation 
6(j) on the gauge field dynamics. Consequently, we find that the maximum value of in 
our framework is 10“^®G and it is not sufficient to explain the Blazar observations. However, 
the above two assumptions can be relaxed, for example, by phenomenologically taking into 
account the effect of charged particles [43] and performing lattice simulations, respectively. 

We have considered the inflaton with the quadratic mass potential and the axial coupling 
in this paper, because it is one of the simplest setups to generate magnetic fields in the 
primordial universe. Unfortunately, however, we have found that the produced magnetic 
fields have a too small correlation length, Aphy(tnow) ^ O.Ipc, and hence its effective strength 
is too weak to explain the blazar observations. Therefore, as a natural extension of this 
work, it would be interesting to consider a scenario in which magnetogenesis due to the 
axial coupling takes place at a lower energy scale. For example, in a small field inflation 
model, inflation ends at a lower energy and hence magnetic fields are produced on a larger 
scale. Alternatively, one can also consider that a spectator (pseudo-)scalar field has an axial 
coupling to the gauge field. When the spectator field begins to oscillate during the oscillation 
phase of the inflaton or (dark) radiation dominated era after inflation, electromagnetic fields 
can be generated in a similar manner to that found in this paper, if the inflaton decays 
into a dark sector without electrically charged particles. Although models of a low energy 
inflation and/or a spectator field contain more parameters than the simple case explored in 
this paper, one can compute the resultant strength of the magnetic fields in the same way, 
once a concrete model is specified. 
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